CONVERGENCE OF NUMERICAL SCHEMES FOR SHORT 
WAVE LONG WAVE INTERACTION EQUATIONS 
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Abstract. We consider the numerical approximation of a system of partial 
differential equations involving a nonlinear Schrodinger equation coupled with 
(^J a hyperbolic conservation law. This system arises in models for the interac- 

tion of short and long waves. Using the compensated compactness method, 
we prove convergence of approximate solutions generated by semi-discrete fi- 
nite volume type methods towards the unique entropy solution of the Cauchy 
£N| problem. Some numerical examples are presented. 

< 

1.1. Interaction equations of short and long waves. The nonlinear interac- 
ts 

tion between short waves and long waves has been studied in a variety of physical 
situations. In [4j, D.J. Benney presents a general theory, deriving nonlinear differ- 
ential systems involving both short and long waves. The short waves u(x, t) are 
described by a nonlinear Schrodinger equation and the long waves v(x,t) satisfy a 
quasilinear wave equation, eventually with a dispersive term. In its most general 
f" — form, the interaction is described by the nonlinear system 

(N 

| idtii + ic\d x u + d xx u = au v + j\u\ u 
\ d t v + c 2 d x v + ndlv + vd x v 2 = I3d x (\u\ 2 ), 
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<N (1-1) 



where c%, c%, a, ft, 7, \i and v are real constants. 

Systems of this type have attracted attention ever since the pioneering works of 



ON 

M. Tsutsumi and S. Hatano [501 [3T], wn0 set c<z = fj, = v = in (1.1 1 and proved 
# ^ the well-posedness of the Cauchy problem. This result was later generalized by 

Bekiranov et al. 3 . More recently, there has been renewed interest in systems of 
5_i the type Here we just refer to [S] for the study of a Schrodinger-KdV system, 

and to 10 for an in-depth study of several generalizations of (1.1), among which 
is the one concerning us in this work. 



Regarding the numerical study of systems of type ( 1.1 ), the only available works 
are, to the authors' knowledge, the papers pQ and [7]. In pQ, the convergence (in 
H 1 x L 2 ) of a semi-discrete finite difference approximation of a version of the system 



( 1.1 ) with C2 = /i = v = is proved. 
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1.1.1. A Schr 5 dinger -conservation law system. Here, following Dias et al. |10j . we 
consider the case where the long waves are modelled by a nonlinear hyperbolic 
conservation law. Thus in [10] the following system is proposed, 



(1.2a) 
(1.2b) 
(1.2c) 



idtU + d xx u = \u\ 2 u + ag(v)u 
d t v + d x f(v)=ad x (g'(v)\u\ 2 ) 
u(x, 0) = u (x), v(x, 0) = v (x), 



where / £ C 2 (M.) and g £ C 3 (K) are real functions such that /(0) = and g' has 
compact support, and a > is a constant. In that work, the authors establish the 



well-posedness of the Cauchy problem ( 1.2a)-( 1.2c I, for an appropriate notion of 



entropy solution, provided a is small enough. 

Comparing this model to the original formulation of Benney, we remark the 
appearance of a general interaction function g(v) in the Schrodinger equation and 
of its derivative g'(v) in the conservation law. This new function, it is argued in 
(10j . gives a more physically realistic coupling between the two equations. The fact 
that g' has compact support allows the authors of [10 to obtain a uniform bound 
for Halloo through a sort of maximum principle. This L°° bound is fundamental in 
the analysis, and it is not clear whether a corresponding well-posedness result can 



be obtained for system ( 1.2a |-( 1.2c) with g(v) = v, and for general / (see, however 
[§] and [2] for some partial well-posedness and blow-up results, respectively). 

Our objective here is to extend the numerical results in [1] to deal with the system 
( |1.2a| -( |l~2c| ). We will prove the convergence of a large class of schemes using the 
compensated compactness method |11[ 118] , aided by energy inequalities which play 
the role of the classical entropy inequalities in the analysis of the method. 



We may summarize some of the main difficulties as follows. The equation (1.2b) 



may be seen as a conservation law with a non-homogeneous flux-function given 
by f(v) — g'(v)\u\ 2 , or as a conservation law with the source term d x (g' (y)\u\ 2 ). 
However, to the best of our knowledge, the available convergence techniques for 
this kind of problem would require a pointwise bound on \u\ 2 (see, for instance, 
[HI E] ) - For problem (1.2a)-(1.2c), such a bound is only obtained a posteriori, and 
only as a consequence of the fact that u £ H 1 (M.) and a Sobolev imbedding. 

Therefore, the crucial bounds in Lemma |3.3| and Proposition |3.4| below must em- 
ploy new techniques which take into account the nonlinear coupling of the equations 



(1.2a ), (1.2b ), as was already observed in ffi HP]. In the present paper we combine 
the techniques of [5J [TU] with those of the compensated compactness method to 
derive new, stronger, compactness estimates which ensure the convergence of our 
numerical schemes. 



1.2. Outline of the paper. In this paper, we establish the convergence of a class 
of numerical schemes to approximate the problem ( 1.2a|— ( 1.2c |. To this end, we 



use a semidiscrete finite volume type scheme as an approximation of the quasilinear 



equation (1.2b) and a semidiscrete finite-difference scheme for the nonlinear Sch- 
rodinger equation, (1.2a). The compensated compactness method allows us to 



prove convergence towards the unique entropy solution, a result which had been 
announced (for the particular case of the Lax-Friedrichs scheme) in pQ. 



The convergence result concerning problem ( 1.2a)-( 1.2c ) (Theorem 2.2 below) 



relies on the compensated compactness method and on the crucial energy estimates 
in Lemma |3.3| and Proposition |3.4| Interestingly, these estimates are finer than 
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the corresponding ones in |10j . and as a consequence our result does not require 
a smallness assumption on the coupling parameter a in (1.2a I, (1.2b I, and which 



we take equal to unity for simplicity. Also, no existence of solution is assumed a 
priori, and thus our convergence proof is also a new existence proof for the problem 



(1.2a|-(1.2cl 



An outline of the paper follows. In Sect ion [2] we define our numerical method and 
state the main convergence result (Theorem |2.2[ ). In Section [3] we prove the main 
stability and energy estimates. In Section[4]we apply the compensated compactness 
method and the estimates of Section [3] to deduce compactness of the approximate 
solutions. We then prove our convergence result. Finally, in Section [5] we present 
some numerical results and discuss some open questions. 



2. The numerical scheme and statement of the main result 



First of all, we recall from [TU] the notion of entropy solution to problem ( 1.2a I 



(1.2c| 



Definition 2.1. Let rj(v) be a convex function (the entropy), and define the entropy 
fluxes <2i,<72 by q[(v) = rf(v)f'(v) and q' 2 (v) = rf(v)g"{v). We say that (u,v) £ 
Lg > c (M x [0, oo)) is an entropy solution to the problem (1.2a)-(1.2c) (with a = 1) if 
for each entropy triplet (77, 91, 92) we have: 

(1) u £ L^ c {[Q,oo)\H x (R)) n C7([0,oo);L 2 (M)), u(Q) = u in L 2 {R), and 

iud t 9 + d x ud x 9 + {\u\ 2 u + g(v)u)9 dxdt = 

tx[0,oo) 

for every 6 € Cg°(R x (0,oo)); 

(2) For every non-negative <f> = 4>(x,t) € C§°(M. 2 ), 



x[0,oo) 



r)(v)d t 4>+( qi (v)-q 2 (v)\u\ 2 )d x <t> 



+ {i 1 '{v)g'{v)~q 2 {v))d x \u\ 2 (t)dxdt+ / r](v o (x))<j)(x,0)dx>0. 
2.1. Finite volume schemes. In this section, we define a class of semidiscrete 



finite volume schemes to approximate ( 1.2a)-( 1.2c I 



Consider a uniform spatial grid (xj)j^ of element size h 



L 3 + l 



xa, where 



h will be the discretization parameter. We seek functions (u h (t) , v h (t)) where for 
each t > u h (t) and v h {t) are sequences u Jt) € C and Vj(t) £ M. 

For the nonlinear Schrodinger equation (1.2a I, we use a standard finite difference 



scheme. The term d xx u is approximated by the usual discrete laplacian, denoted 

by 

1 

)] = U h u 3 - J^2 

We will also need the difference quotient 



(Dlu h )j = D 2 hUj = ~^{u j+1 - 2uj + Uj-x). 



D + Uj 



For the second equation, (1.2b), we will consider a finite volume scheme, as 
follows. First, let us rewrite ( 1.2b[ ) in the form 

(2.1) 



d t v + d x (f{v)-g'(v)\u\ 2 ) =0. 
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Now for each j € Z consider monotone, conservative, Lipschitz continuous numeri- 
cal flux-functions 

f 3 ,+ (vi,v 2 ), fj-(vi,v 2 ) 

consistent with the function /, that is, fj,+ , fj— verify for each j € Z, V\, v 2 , v £ R, 
the conditions 

(2.2a) fj,+ (vi,v 2 ) = -fj+i t -(v2,v\) (Conservation) 

(2.2b) f jt ±(v,v) = ±f(v) (Consistency) 

(2.2c) difj,± > 0, 9 2 fj,± < (Monotonicity), 

where d\ and d 2 denote the derivative with respect to the first and second argu- 
ments, respectively. 

Next, consider another family Gj !+ (vi,v 2 ), Gj,-(vi,v 2 ) of numerical flux-functions 
verifying the conditions in (2.2) but with —g'(v) instead of f(v). Thus, 

(2.3a) Gj t+ (vi,v 2 ) = -Gj+i t -(v2,vi) 

(2.3b) G ji± (v,v)=T9'(v) 
(2.3c) 9i<5 ii± > 0, d 2 G jt ± < 0. 

Further, we consider a family of Lipschitz continuous functions 

Gj,+ (vi,v 2 ,ai,a 2 ), 

satisfying the following conditions: 

(2.4a) G jt+ (vi,v 2 ,ai,a 2 ) = -G j+i _(v 2 ,v i ,a 2 ,a i ) 

(2.4b) G jt ±(v,v,ai,a 2 ) = Tg'(v)-(ai + a 2 ) 

(2.4c) Gj t ±(v 1 ,v 2 ,a,a) = aG jt ±(vx,v 2 ) 

(2.4d) dxG j>± > 0, d 2 G jt ± < 0. 

Thus, if for instance a — \u\ 2 for some complex number u, we find 

G j ,+ (v > v > \u\\\u\ 2 ) = -g'(v)\u\ 2 , 

so that this term gives a consistent approximation of the term —g'(v)\u\ 2 . 
Finally, define Hj ± by 

(2.5) Hj,±(vi,v 2 ,ai,a 2 ) = fj,±(vi,v 2 ) + G jt ±(v 1 ,v 2 ,ai,a 2 ). 

In short, Hj± are Lipschitz continuous, monotone, conservative, numerical flux 
functions, consistent with f(v) — g'(v)\u\ 2 in the sense that 

(2 - 6) Hj t ±(v, v, \u\ 2 , \u\ 2 ) = ±f(v) + Gj,±(v, v, \u\ 2 , \u\ 2 ) = ±f(v) T g'(v)\u\ 2 . 

For background on finite volume schemes, see [T2l 114) . 

We now approximate the equation ( 2.1 1 by a semi-discrete finite volume method. 
We give some initial data uoji^oj on < = 0; for instance, uoj — ^ J X J+1 uo(x)dx 
(and similarly for w ), so that the piecewise constant functions defined from uqj, 
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voj converge strongly in Lj oc (R.) to uo,vo. We then solve the infinite system of 
ODEs 



(2.7a) idtiij + D h uj = \uj\ Uj + g(vj)uj 
1 
h 

2.2. Examples 



(2.7b) d t Vj + \{H Jt+ (v h v j+u \u 3 \ 2 , \u 3+1 \ 2 ) + H jt _(Uj,Uj_i, \u 3 \ 2 , |^j-i| 2 )) = 0. 



2.2.1. A Lax-Friedrichs scheme. Perhaps the simplest example of a numerical flux 
satisfying the above conditions is given by the following Lax-Friedrichs type scheme. 
Let A, 7 be some constants and set 

(2.8) 

fj,±(vi,v 2 ) = f±(vi,v 2 ) = ±-(/(wi) + f(v 2 )) T ^(w 2 - vx), 
G 3 ,±(vi,v 2 ,ai,a 2 ) = G±(vi, v 2 , at, a 2 ) 

= -g'{vi)ai - g'{v 2 )a 2 ) =F ^{v 2 - ui)-(ffli +a 2 ). 

It is easy to check that, under an appropriate CFL-type condition on A and 7, the 
numerical fluxes f± and G± verify the assumptions of the previous section. 

2.2.2. A Godunov scheme. Here we consider the following straightforward general- 
ization of the classical Godunov scheme. Setting 



Hj,± = H±(v 1 ,v 2 ,ai,a 2 ) 



min ±(f(s)-g'{s)ha 1 +a 2 )), v x < v 2 , 

Vi<S<V2 Z 

max ±(/(s) - g'(s)-(ai + a 2 )), v 2 <Vi, 

V-2<S<Vi Z 



one can then verify that the consistency, conservation and monotonicity conditions 
above are met. 

2.3. Convergence result. We are now ready to state the main result of this paper, 
whose proof will be the object of Section [4j 

In what follows, we suppose that / and g verify a standard non-degeneracy 
condition: 

Vk>0, the set {s : f"(s) - Kg"'{s) ^ 0} is dense in R. 
This condition is needed in order to apply the compensated compactness method. 



Theorem 2.2. Let (u h ,v tl ) be defined by a semidiscrete scheme (12/7]), verifying all 



the assumptions (2.2)-rt2.4l. Then there exist functions 



u6£7([0,c»);ff x (R)), v e L°°{R x [0, 00)), 



solutio ns o f the Cauchy problem ( |1.2a |-( 'l.2c[) (with a = 1 ) in the sense of Defi- 



nition 
[0,oo)J 



2.1 



such that, up to a subsequence, (u h ,v h ) converge to (u,v) in £j oc (M x 
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3. Uniform bounds and energy estimates 

In this section we prove the crucial uniform estimates of Propositions |3.2[ |3.4 
and Lemma |3.3| We begin with a key property of finite volume schemes. 

Let be a numerical flux-function consistent with a function /, let V\, v 2 £ K 
and let r\ be a smooth convex function. Define the viscosity as the quantity 



(3.1) / r ] "(v)(f(v)-f j , + (v 1 ,v 2 ))dv. 

J vi 

The following property, although classical, is very important in the analysis of the 
scheme (cf. [HIE]). ^ * s an easy conse q uence of the monotonicity and consistency 
properties of the flux functions, and so we omit its proof. 

Lemma 3.1. For each Vi,v 2 £ K and each smooth convex function, the viscosity 



(3.1) is everywhere non-negative. In particular, we have 

V / '(v)(f(v)-f J . :+ (v 1 ,v 2 ))dv>0, 



V2 

n"{v){-g'(v) - G j>+ (vi,v 2 ))dv > 

for every v\,v 2 £ K.. 

In what follows, if (say) v h is a sequence, we will not distinguish between the 
discrete norm of v h in l p (Z) spaces its and continuous norm Hi^Hp defined as the 
usual LP norm of the piecewise constant function v (x) — Vj if x £ [xj, Xj+i]. Thus 
if p £ [l,oo), 

ik' i iip-(E^i p ) 1/P ' 

with the usual modification if p = oo. 

We now prove an L°° bound for v h which is essential in all our analysis, as well 
as the conservation of the L? norm of u . 



Proposition 3.2. Let (u ,v ) be defined by the finite volume method (2.7), for 
some initial data uqj G I (Z), VQj £ I (Z). Then, for each h > 0, there is a unique 



global solution of (2.7). Moreover, for some M = M(vo,g) > independent of h 
(where g is the coupling function in ( |1.2a[ ) ), we have 

(3-2) j t \\u h {t)\\l=0 
and 

(3.3) II^WIloo < M. 



Proof. First, the local in time existence of solution of (2.7) is ensured by a standard 
fixed point argument in the discrete space Z 2 (Z). The argument is very similar to 
the one in [TU], and so we omit the details. 



Next, we sketch the proof of estimate (3.2), which consists of a straightforward 



calculation. From the scheme (|2.7a|), multiply by huj, sum in j £ Z and sum by 



parts the term with D^Uj to obtain that 



idt^hluA 2 £ 



"3\ 
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which immediately gives (3.2) 



Since, for fixed h > 0, one has P(Z) C Z°°(Z), the local solution will be global in 



time if we prove the bound (3.3) 



Following [TU], consider M' such that supp g' C [-M',M'\. Then, setting M = 
max{M', ||i>' l (0)|| oo }, we will see that \\v \\oo < M. First, we consider a perturbed 
problem which coincides with (2.7), except that the term — esgnv h is added to the 



right-hand side of the second equation, for some small e > 0. Let (u h,e ,v h,e ) be the 
solution of this problem, whose existence is obtained in the same way as for the 
unperturbed problem. 

In what follows, let us consider only the terms involving the coupling function g' , 
since the treatment of the terms involving the flux / is standard and can be found 
for instance in [14] . Thus, suppose that there is a first j 6 Z and a first instant 
t > such that Vj(t) = M (the case Vj(t) = — M is similar). Then the assumptions 
on the coupling function g, namely the compact support of g' , imply that the term 



with Gj.± appearing in the definition of the scheme (2.7b) may be written as 



1 



- -(G ii+ (^ +1 , K-l > K+il 2 ) + G,,-(vh v U> Kt> K--il 2 )) 



;G,, + (^ +1 , \uf, \u) +1 \ 2 ) + G jt+ (vl,vl, |uj| a , |^ +1 | 2 )) 
(Gj,-( v P v j-i> K-l 2 , K-il 2 ) + Gj- ( v h v h K-l 2 ' l^i-il 2 )), 



due to g'(vj) — and the consistency conditions in ( |2.3[ ),(2.4). Now divide and 
multiply the first sum by — (which is non-positive since v ' e attains its 
maximum value at j), and the second sum by Vj_i~ v^. The monotonicity condition 



in ( 2.4 ) then ensures that the resulting terms have the right sign so that (along with 



the treatment of the terms involving fj +, which we omit 



dtv^t) < -e < 0, 

which is absurd, since we assumed that v ' e attained a maximum value at j, t. Thus, 
\v h ' e \ < M. 

It is now necessary to take the limit e — > (with h fixed) to deduce the de- 
sired estimate (3.3). We obtain from the scheme ( 2.7b[ ) using the regularity of the 



numerical flux functions H 



(3.4) sup \v e j{t) - Vj(t)\ < C{h) I sup \Vi 
j£Z Jo j& 



sup \u 



e 1 2 



ds + et, 



where C(h) is some (unbounded) function of h which may change from line to line. 
Note that the estimate (3.2) remains valid with e > with exactly the same 



proof, that is, we have ||u' l ' e (t)||2 = ||w' I (£)||2 = C. From I 2 c l°° we get immediately 
lK' e ||oo < ft- 1/2 ||tt ft ' e ||2- Thus, from J3~4|, 



sup I \u) | 2 

j 



< 2 sup |Uj||uj - 

3 

< 2/i- 1 / 2 ||m' 1 || 2 su P | 



Ua — Uj 



< C(h) sup \u) 
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and using the scheme (2.7a) gives in turn 



sup \uj —Uj\ < C{h) \ sup \uj — Uj\ + sup \vj — Vj\ ds. 

jez Jo jez, jez 



Gronwall's lemma yields 



sup \Uj — uj\ < C(h, t) sup \Vj 



which we substitute in (3.4) to get (after applying Gronwall again) 

sup|^(t) - Vj(t)\ < C(h,t)e 



which goes to zero with e, for h and t fixed. Therefore we deduce (3.3). This 
completes the proof of Proposition |3.2| □ 



We now prove an energy inequality for the second equation of (2.7), involving 



the viscosity (3.1), which will be essential to establish Proposition 3.4 



Lemma 3.3. Let (u h ,v h ) be defined by the scheme (2.7) and let i] be a smooth 
convex function. Then, the estimate 

Hv h (t))\\i+ r + \"(v)(f(v)-f j , + (vj,v j+ i))dv 



(3.5) 



<c + c 

holds for some constant c > independent of h. 



v"(v)(-g'{v) - Gj.+ivj^j+^dvds 

f \\D + u h {s)\\ 2 ds. 
Jo 



Proof. The scheme (2.7) reads 

+ j^(Gj,+{'"j,v j+1 , \uj\ 2 , \u j+1 \ 2 ) + G jt -(vj,Vj-u M 2 , \u 3 -i\ 2 )) = 0. 

Now multiply by hrj'(vj) and sum over j € Z to obtain, with obvious notation, 

d t ^2hri(v j ) + ^2ri'(v j )(f j>+ (v j ,v j+1 ) + f :i ^(v j ,v j _ 1 )) 
jez jez 

+ ^2 r )'(vj)(Gj, + {v j ,v j+ i, |^| 2 , \u j+1 \ 2 ) + Gj-(vj,Vj-i, \uj\ 2 , |uj-i| 2 )) 

ee d t E hr]{Vj) + A + B = Q. 

jez 

First, summation by parts gives 

A = -22(v'(v j+ i)-r)'(v :i ))fj i+ (y j ,v j+1 ) = -22 "n"{v)f jt+ {vj,v j+1 )dv. 
jez jez v i 

Next, let qi(v) be such that frf = q[. We have 

o = E , ?i(^+i)-9i(^) = E / iWW* = -E v "(v)f(v)dv, 

jez jez ^ jez v i 
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and therefore 

(3.6) A = ~z2 V"{v){f jl+ (vj,v j+1 ) -/(«)) dv. 

jeZ Jv i 

Now consider the term B. First, we find 

B = - ^2(r)'(vj+i) - r 1 '{v 3 ))G j .+{v J1 v J+ll \u 3 \ 2 , \u j+1 \ 2 ) 

jez 

rvj + i 

= ~Z2 v"i^)Gj,+(vj,v j+1 , Kl 2 > \u j+ i\ 2 )dv. 
jez ■' v 3 

Next, setting q' 2 — rfg" gives 

= E 92 (Uj+i ) | %+i | 2 - 92 ) |«j 1 2 

= J2(l2(vj+i) - q 2 (vj))\uj\ 2 + q 2 (v j+ i)(\u j+1 \ 2 - \uj\ 2 ) 
jez 

= El^! 2 9"(v)r,'{v)dv + q 2 (v j+1 )(\u j+l \ 2 - \ Uj \ 2 ), 

jez Jv i 

and thus, integrating by parts, 
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jez uv i jei 
J2q2(v J+ i){\u J+ i\ 2 - \uj\ 2 ) 



= - EN' / 9'(vW'(v)dv-(q 2 (v j+1 )-g'(v j+1 W(v j+1 ))(\u j+1 \ 2 - \ Uj \ 2 ). 
jez J ^ 

This gives 

B = -E r ,+ V , («)(G f i,+(t; J ,» J+1 ,|u i | a ,|« j+1 | 2 ) + |u i |V(i;))d« 

jeZ Jl, J 

+ E M^j+i) -5'(^'+i)?7'(^+i))(l«i+i| 2 - M 2 ) 
3'ez 

= _ E l u J'| 2 / 7 /'( u )( g j\+(^ u j+i) + .9») d« 

jez 

+E (feOj+i) -5'(^+i)??'(vj+i))(i«j+ii 2 - n 2 ) 

jez 

+ E / ^'( w ) du (kl 2G .7\+( v i> w j+i) -Gj, + (Wj,Wj+i,|%| 2 ,K- + i| 2 )J. 
jez- 7 ^ 

The first sum in the right-hand side of the previous equality is the sum appearing 
in (3.5 1. The other two sums are bounded by 

c£|kf-h+il 2 | 

jez 
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by appealing to \uj\ 2 G jt+ (yj, v j+ i) = G jt+ (vj, v j+ i, |u,-| 2 , \uj\ 2 ) (cf . (pl| ), the 
regularity of the numerical flux functions Gj :+ , the L°° bound on Vj, (3.3), and the 
smoothness of g and r\. Finally, in view of the bound on the L 2 norm of u , (3.2 1, 
we find 



El 



\uj\ 2 - |uj+i| 2 | 



< 



(3.7) 



Thus 



\u j+ i - UjWuj^+uj] 

<2(^iK +1 -u J f^/i|u j 
= 2\\D + u%\\v h \\ 2 = c\\D + u%. 



1/2 



B = E 



i£2 



with < c||Z? + u ||2- This estimate and (3.6) give (3.5), after integration on 
(0, t). This completes the proof of Lemma 3.3 



□ 



The following is the crucial a priori estimate in our analysis. 



Proposition 3.4. Let (u h ,v h ) be defined by the scheme (2.7). Then, there exist 
non-negative functions a(t), b(t) continuous on [0, oo) such that 



(3.8) 
(3.9) 



||t> h (*)|| 2 <&(*), 
\\D+u h {t)\\ 2 < a(t). 



Remark 3.5. In view of the estimate (3.9), the viscosity estimate (3.5) implies 

ft 



(3.10) 



o 



E 



Vj + l 



v"{v){j{v) - f j:+ (vj,v j+1 ))dv ds 



Vj+l 



rf'(v)(-g'{v) - Gj t+ (vj,v j+ i))dvds < c(t) 
for any smooth convex function r\ and some locally bounded function c(t). 



Remark 3.6. In the case of the Lax-Friedrichs scheme (2.8), the estimate (3.10) 
actually implies an explicit bound on the "quadratic total variation" of v , as was 
observed in jT]. More precisely, we have for some k > 



(3.11) 



i« h (*)in+* Te^ 



+ \uj\ )(vj + i — vj) 2 ds < c(t). 



This can be seen by observing that for the Lax-Friedrichs scheme, one has (for 
i>i < v < v 2 \ the other case is similar) 

f( v ) - U(vi,v2) f+(v,v)-f + (vt,v 2 ) 



V 2 -Vl V 2 ~ V\ 

f+(vi,v) - f+(vi,v 2 ) _ 
v 2 - Vl 

f + {vi,V 2 ) - f+{vi,V 2 ) 
V 2 - Vl 



> 



> 



k = k>0, 
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for some k > 0. This follows from the monotonicity of the scheme, and k appears 
due to the particular form of the Lax-Friedrichs scheme and a CFL condition 
A|/'| < 1. Thus (see plo| , with if = 1) 

r>->2 

(3.12) 



f(v)-i + (v 2 ,vi)dv> / k(v 2 - vi) dv = k(v 2 - «i) 



This, together with a similar estimate for the numerical flux G + in (2.8) gives 
(3.11). In the general case, however, the estimate (3.12) is only valid for k = 
(for example, the Godunov scheme in the neighborhood of a sonic point fails to 
verify (3.12) with k > 0). Thus, in general there is no bound on the quadratic total 
variation as in ( 3.11[ ), and the appropriate estimate is (3.10). 

Proof of Proposition^^ First, of all, recall the Gagliardo-Niremberg-Sobolev in- 
equalities which we will use throughout. If <j> € Z 2 (Z), then 



(3.13) 
(3.14) 



1/2 1 



1 1/2 
2 



U\\oo<C\\<t>\ 

U\U<c\\<j>t 2 /A \\D + ct>\\y\ 



Now take (2.7a), multiply by hd t Uj and sum over j £ Z to obtain, using straight- 
forward computations very similar to the ones detailed in [101 Lemma 2.2], 



(3 



15) l\\D + u% + \\\u% + lj2 h 9( v ^\ 2=c + / E^'(«i)^l % | 2 d S . 

JO jgZ 



Let us now consider the term on the right-hand side of (3.15). Using the second 
equation of the scheme (2.7b), and the conservation properties of the numerical flux 
functions, wc find (omitting the arguments of the numerical flux functions) 



^hg'(vj)d t v 3 



Kf (/*+ + h- + G jt+ + G jt -) 



(3.16) 



E(5'fe+i)K+i| 2 -9'fe)Kf) (f jt+ +G j>+ ). 



First, we have 

E {9'{v ]+ i)\u J+ i\ 2 - g'{v j )\u j \ 2 )f j .+ = \uj\ 2 {g'(v j+ i) - g'{vj))fj t + 



(3.17) 



.ie- 



.ie- 



Efl'(^i)(i' 



A 1 + A 2 . 



We will write A 1 in terms of the viscosity in order to use the estimate (3.10). 
^i=EM 2 P +1 9"(v)f 3 , + dv 

= EN 2 9"(v)(f j!+ -/(«)) d^ + EKf 9"(v)f(v)dv. 

It is an easy consequence of the monotonicity and consistency properties of the 
scheme that for s between Vj and Vj+i, the quantity fj : +{vj, Vj+i) — f(s) does not 



12 



P. AMORIM AND M. FIGUEIRA 



change sign (cf. Lemma 3.1). Therefore, we may write for some intermediate value 
c (observing that from Lemma 3.1 the integral below has the right sign) 



g"(v)(U + - f(v))dv = \g"(c)\ / (/(«) - h + )dv 



Moreover, defining a function 3" by = g"(v)f(v), we find 

E N 2 9"(v)f(v) dw = E l%| 2 (^i+l) - Hvj)) dv 

= - E (i^+ii 2 ~ n 2 ) dv i 



and thus, by the uniform bound for v h (3.3) and the Sobolev inequality (3.13), 

/Vj + l 
(/(«) - /,,+) du + E sup |5] I 1 2 - M 2 | 

< c||£> + u*|| a E (/(«) - d« + c||D + u h || a . 
We may now apply Lemma 3.3 with rj(v) — v 2 /2 to get after integration on (0, t) 

Ai(s)ds < sup||L> + u' l (s)|j 2 (c + c / ||L> + M' I (s)|| 2 ds) 
(0,t) v Jo ' 

+ c / ||D + u h («)|| 2 <k. 



The term A 2 in (3.17) is estimated more simply using the uniform bound for v h 
(3.3) to give after integration on (0,t) 

A 2 (s)ds <c [ \\D + u h (s)\\ 2 ds. 



Thus, (3.17) becomes 



(3.18) 



< 



(CO 



fj,+ds 








J \\D + u%dsJ 


+ c f \\D + u h \\ 2 ds 




Jq 



The remaining term in (3.16) gives 



(3.19) 



E (s>, + i)K+i| 2 " ff'(«j)l« J f)G ij+ = E M V(«i+i) - fl / («i))G i>+ 
jez jez 

+ Eff'^+i)(l^+i| 2 - l^'| 2 ) G i.+ = : Bi+B 2 . 



We have 



(3.20) 



B 1 =J2\u j \ A (g'(v j+l )-g'(v j ))G j ,. 



EK-| 2 (9'(^+i)-5'(«i))(%+ ~ \Uj\ 2 G jt+ ) =: Bn + B 12 . 
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Again, we use estimate (3.101, and so we need the viscosity to appear. We find 

rvj+i 



g"{v){G j ,+ +g'{v))dv-Y,M i / g"(v)g'(v)dv 
< ||^||LE su Plff"IKI 2 / (-G j<+ -g'(v))dv 

+Ei(-9'(^+i)) 2 lK+ii' 



14 luil 4 ! 



which gives, using Lemma 3.3 and the Sobolev inequality (3.131 

Buds < c sup \\D+u h \\ 2 I 22\uj\ 2 I (~Gj,+ - g'{v))dvds 



(o,t) 



+i| 2 - |uj| 2 Ms 



< sup \\D + u h \\ 2 (c + c f ||U + u & ||ad8 
(0,t) v Jo 



The term f*Bi 2 ds from (3.201 yields using (3.3) and ( 3.13[ ) (recall the relation 
between and \v,j\ 2 Gj. + in \2A\) 



[ B 12 ds< ( |mj| 2 LipG||w J+ i| 2 - \u 3 \ 2 \ds < c sap \\D + u% [ \\D + u h \\ 2 ds 

JO Jo - p y (0,t) Jo 



As for B 2 in (3.19), wc find immediately using the uniform bound for v h (3.3) 



[ B 2 ds<c [ \\D + u h (s)\\ 2 ds. 
Jo Jo 



Thus, (3.16) becomes 
ft 



(3.21) 



^hg'iv^dtvjlujfds 

jez 

<sup||D + « fc || 2 fc + c / \\D + u%ds) + c f \\D + u%ds, 

(0,t) V I" ' In 



and, after observing that estimates rt3.2h and (3.3) imply 



we conclude from (3.15), (3.18) and (3.21) that 



\\\D + u h (t)\\l + \\\u h (t)\\i < c + c jf 1 \\D + u h \\ 2 ds 

+ sup \\D+u h \\ 2 (c + c [ 
(o,t) v Jo 



|D + tt ft || 2 ds 
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In particular, setting /3(f) := sup( t ) (l + \\D + u h (t)\\f) this easily implies 

/3(f) < c + c f t3{s) 1/2 ds 
Jo 

+ /3(f) 1 / 2 (c + c J PisY^ds}. 
The estimate (3.9) follows by Gronwall's lemma applied to the function /3(f) 1 / 2 . To 



prove (3.8 1, simply apply (3.9 ) to the estimate (3.5 1 in Lemma 3.3 with tj(v) — v /2 



This completes the proof of Proposition |3.2 



□ 



4. Proof of convergence 



In this section we prove Theorem |2.2| relying on the compensated compactness 
method [TTJ [TS] . According to this method, the strong compactness of a sequence 
of approximate solutions iy h ) is a consequence of the following property: 

(4.1) d tV (v h ) + d x { qi {v h ) - \u h \ 2 q2{v h )) € { compact of W^' 2 }, 

where, we recall, 77(f) is a convex function (the entropy), and the entropy fluxes 
<Zi.2 verify q[(v) = T]'(v)f'(v) and q' 2 {v) — T)'{v)g"{v). In practice, one may use the 
following well-known result to establish (|4.1[): Ifl<<7<2<r<oo, then 



(4.2) { compact of WZt A } n { bounded in WZ}' r } C { compact of WZ}' 2 } 



Lemma 4.1. Let (u , v ) be defined by the semidiscrete approximation (2.7 1. Then, 



the compactness property in (4.1) is valid. 



Proof. Let 4> be continuous and compactly supported on K x (0, 00). Let us fix some 
notations. In what follows, <f>,j = cf>(xj,t), Ij — [xj, Xj+x] and = ^ Jj <f>(y,t)dy. 
Let J = J(h) e N, T > be such that supp</> C [x-j,xj] x [0,T]. Also, we 
sometimes write, say, r/L := T)'(vj), and likewise for the other functions. We view 
v and u h as piecewise constant functions defined as v (x) = Vj if x £ Ij . We have 



from (4.1 ) 



(4.3) - (d t n{v h ) + d x (qi(v h ) - q2(v h )\u h \ 2 ) 1 4>) v , xV 



hy^fyr)' (vj)d t Vjdt - 

n j& JO 



E 



4>j+i(qi{vj+i) - qi{vj))dt 



E 



(f>j+i(q2{vi+i)\u 



q2(V])\u 3 \ 2 )dt, 



where we have just differentiated piecewise constant functions in the sense of dis- 
tributions. From the definition of v h in the scheme (2.7), we find, with obvious 
notation, 



fijri'iv^dtVjdt = 



4>jV'{v ] )(f h+ {v ]1 v ]+1 ) + fj-{vi,Vj-x))dt 







^2 ( />jV'(v j )(G jj+ {v j ,v j+1 , \uj\ 2 , \u j+1 \ 2 ) + G jt -{vj,Vj- U \uj\ 2 , |uj_i| 2 ))df 



=:A + C. 
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Thus, if wc define (see (gjjj)) 

pOO 

B :=- y2<f> j+1 (q 1 (v j+1 ) - qi(vj))dt 
Ja jew, 

pOO 

D := / ^^j+ii^ivj+i^Uj+il 2 ~ q 2 (v J )\u j \ 2 ) dt, 

JO £ r* 



we get 
(4.4) 



- (ftx, t)d tV (v h ) + d>(x, t)d x ( qi (v h ) - q 2 (v h )\u h \ 2 )) v , x v 
= A + B + C + D. 

We must therefore estimate A + B + C + D so that the compactness property in 
( |4.1[ ) is valid. 

Estimate of A + B. We begin with (omitting the arguments of the numerical flux 
functions whenever possible, so that fj t + = fj t +(Vj,Vj+\), fj - = fj,-(vj,Vj-i)) 

poo poo 

i / E + h-) dt + / - & + h-) '" 

=: Ai + ^la- 



First, using the conservation property in (2.2) and summing by parts, we find 



/>OG 

POO POO 

= ~ I - 'PjWjfj^ dt- E ^+1(^+1 - Vj)fj,+ dt 

Jo ^1 Jo „.,-„ 



and so 



pOO poo f v j + l 

Ai = - £> 3 -+i - ^Hh+ dt- E fc+l / 7?" , 

JO JO J-U,- 



=: A n +A X2 . 
Next, we have 

B = - E^+i / V'(v)f(v)dvdt 

JO j& Jvj 

= E^+i / v"W(v)dvdt- Y,^+M+if3+i-^fj)dt 

JO j& Jvj Jo j& 

='■ B\ + B 2 . 
We now write 

(4.5) A + B = T 1 +T 2 +T 3 , 

with 

T 1 =A 12 +B 1) T 2 =A n +B 2 T 3 = A 2 
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and we estimate each term in turn. We have immediately 



|Til< 



E 



V 3 + l 



V"(v)(f(v)-h+)dvdt, 



since by Lemma |3.1|the quantity inside the integral is non-negative. Thus, by the 
viscosity estimate (3.10) we see that T\ is bounded in the space of measures (on the 
support of 0), which is compactly embedded in W / ~ 1,9 (supp (j>) for q £ [1,2). This 
proves the compactness property (4.1 1 for T\. Note also that from Lemma 3.1 we 
have 



(4.6) 



Ti > if <t> > 0, 



which will be of use later. 

We now state a technical result [12l Lemma 4.5], which is essential in our con- 
vergence analysis. Adapted to our setting, it states that if b : K — > K is a monotone 
function, then for some C > 0, we have the bound 



(b(p)-b(q)) 2 <C / b(s)-b(p)ds 



Since f jt+ (p,p) = f(p), this gives 



(f(p) - f j!+ (p, q)) 2 < C / U+(P; s) - h+(P, l) ds 



But the monotonicity of the scheme implies that the integral inside the absolute 
value is actually non-negative, and so 

2 f q 

(lip) - fj,+(p, q)) <C fj,+(p, s) - fj,+(p, q) ds 

(4-7) % 

<C / f(s)-f j ,+ (p,q)ds, 

again by monotonicity. Similarly, we find 

(4.8) (/(g) - .f 3 AP,q)f < C I" f(s) - f j ,+ (p,q)ds. 

Jp 

With an entirely similar proof, we get also for £ — p,q 
(4.9) 



{-g'{i)-G^ + {p,q)Y <C / -g'(s)-G j , + (p,q)ds 



With this result in hand, we now treat the term T2, using again the viscosity 
estimate (3.10). First, some straightforward rearranging gives 

(4.10) T 2 = - / E ^ (^"+1 ~ ^*)(/j,+ - fj) dt - 



Suppose now that <j> is Holder continuous with exponent a £ (1/2,1) and recall 
that supp0 is compact, so that in particular the sum in j takes place over a finite 
set J. We find, using the uniform bound for v h (3.3), the estimate (14. 71 and the 
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property £V £J h = EjeJ Iij 1 dx = Ix-j ldt <C\ supp<£(t, •)!, 

< ch a u\\ 0<a f ( h) 1/2 ( £(4+ - / 3 -) a ) V2 ^ 1/2 * 

<C/l a - 1 / 2 ||^||o.a / (V/ f(v)-f ji+ dv) dt. 



Therefore we may use the viscosity estimate (3.10) to get 
(4.H) \T 2 \<Ch a - l ' 2 \\<j>\\^ a a{T), 

where a(t) is a bounded function. Since W l q C C°' Q , with compact embedding, 
for q' > 2/(1 — a) > 4 (and thus q E (1,4/3)), we see that this term is compact in 
W£« for q e (1,4/3) C (1,2). 

The term T3 is treated in the same way, by putting 

+ fj,-( v j> v j-i) = fj,+ ( v j, v j+i) - f(Vj) - fj-^+ivj-^Vj) + f(vj), 



using (4.8), and observing that \<f>j — <pj\ < /i a ||0||o jC c Thus, the terms A + B in 



(4.4) have the desired compactness property. 



The terms C + D. We now turn to the remaining terms in (4.4). These are non- 
homogenous terms involving the solution u h to the Schrodinger equation, and thus 



present additional difficulties. Recall that C is defined after (4.3). First, we must 
put 4>j in place of 4>j . We have 



C 



poo poo 

JO ^71 JO ■r-v 



far/AG^ + G jt -)dt 



=:C 1 +C 2 . 



Now, from the conservation property of the numerical flux in (2.4) 



poo 

Cl= I X'V/'( ( *'V- :: i- N 2 , \u 3 +i\ 2 ) - Gj-^+iVj-x^j, luj-il 2 , \uj\ 2 )) dt 

JO A rt 



l«ilM«3+l| 2 )* 



pOO 

p OO p OO 

JO -^r, JO 



and so 



poo rcc f v j+l 

ci = - - ^)^g j ,+ - / X! / ??» G ^+ 



=: Cn + C12. 
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On the other hand, 

/■OO 

= / E'AmK+i^M^+i) -q2{vj))dt 

pOO 

+ / y^'+l*2(Uj)(K-+l| 2 - |Mj| 2 ) dt. 

- 70 jez 



But since 92 = J?'5"> we nn d 



= 2^^+il"j+i| 2 / v'(v)g"(v)dv 
= - H &+iK+i| 2 / V"(v)g'(v) dv + Y^ <f>j+i\uj+i\ 2 (Vj+i9j+i - Vjd'i, 



and so 



d = - E^+il^+il 2 / v"(v)g'(v)dvdt 

Jo j& Jvj 

pOO 



/>oo 

+ / y^j+i92( v j)(l«j+i| 2 - |«j| 2 ) d* 



j'ez 

=: £>i +L> 2 +D 3 . 

We now write 

(4.12) C + /J = 5 1 +5 2 + i? 3 + C 2 , 

with 

Si := D\ + C12, 52 := D2 + Cn 
and still (adding and subtracting) 

fVj + l 



(4.13) 



S 1 = - 22<j> j+1 \u j+1 \ 2 T]"(v)(g'(v) + G h+ (v 3 ,v 3+1 ))dvdt 

Jo i6Z J«j 

f 00 \ f vj+1 
+ / z2^+ 1 n"{v){\u j +i\ 2 G jt+ {v j ,v j+1 )-G jt+ )dvdt 
Jo eZ Jt> 3 



j'ez 

=: Sn + Si2- 

Thus we have succeeded in bringing out the viscosity of the scheme in Sn. We now 
bound these two terms. We have 



|Sn| < Halloo / EK+il 2 / rj"^ - g'(v) - G j!+ ( Vj ,v j+1 )) dvdt 

<cu\\c 
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in view of the viscosity bound (3.10) (the bound (3.10) actually involves 1 2 and 
not |u J+ i| 2 , but an examination of the proof quickly gives the same result with 
|wj+i| 2 )- Thus, Sn is bounded in the space of measures, which as we have seen in 
the treatment of the term T\ above, is sufficient for our purposes. Additionally, we 
find from Lemma 13. II 

(4.14) Sn > if > 0, 

which will be useful later. 
Next, we have 

fT 

Sl2 < C\ 



E 



\u j+ i\ 2 -\ Uj \ 2 \dt<CUl 



\D + u 



h \dt, 



by using the uniform L°° bound on v , the bound (3.7) and the smoothness of Gj^+ 
(recall that G Ji+ (a, b, u, u) — uGj 7+ (a, b)). In view of the bound on ||D + u' l ||2, (3.9), 
this term is also bounded in the space of measures. 

Consider now the term 52 = D^ + Cn in (4.12). Some easy rearranging, addition 
and subtraction give 



(4.15) 



oMA 2 &jA - Gj,+) - i+1 r^(|« i+1 | 2 - \u : \ 2 )dt. 



Now the first of these three terms is treated in the same way as the term (4.10) 
above, using the viscosity estimate (3.10), the Sobolev inequality (3.13), and the 
technical estimate (4.9|, giving 



<h a Uh, a \W(v h )\ 



|tt h (t)|U2Juil|0£ + G il+ |dt 

3 6./ 



< Ch 



a u\\ , a / T >Aw(x;^) 1/2 (Ei tt ji a (5J + <5 i.+) 2 ) 1/2/1-1/2 dt 

Jo jeJ jeJ 



<Ch a - l '*\\<t>\\ , a 



< Ch a - x ' 2 U\\ Q , a A{T) 



^w(Em 



jeJ 



-g'(v) - Gj^+dv 



1/2 



dt 



for some locally bounded function A(t). This gives compactness in WZ ' q for q e 
(1, 4/3) C (1, 2). The second and third terms may be treated exactly as S12 above, 
and are thus bounded in the space of measures and so have the desired compactness. 
Note also that the term D 3 has precisely the same form and so yields to the same 
analysis. 

It only remains to estimate the term 



(4.16) 



Co = 



which we only sketch, since no new techniques are necessary: just add and subtract 



appropriately to obtain terms with 



\uiVWi+Gj,+) and 



^■)( G j,+ - 
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\uj\ 2 Gj y+ ), which fall into the cases treated before. This completes the proof of 
Lemma 14.11 □ 



Lemma 4.2. Let (u ,v ) be defined by the scheme (2.7). Then, there exist func- 
tions 

u e L^ c ([0,oo);H l (M.)), v e L°°(R x [0,oo)) 
such that (for a subsequence at least) 



and 



i--tv u h -> u in L\ oc {R x [0,oo)) 



D + u h ^d x u in Z,°°([0,oo);Z/(lR)). 



Proof. We have proved the compactness property (4.1). Thus, according to the 
compensated compactness method [18], there is a function v e L 1 1 oc (M x [0, oo)) 
(and thus in x [0,oo)), by the uniform bound (3.3)) such that v h — > v in 

Lj, c (Rx[0,oo)). 

Let us now analyze the compactness of u h . For this, it will be useful to define 
the piecewise interpolators P\ and Pq. The piecewise linear interpolator P\u h is 
the unique continuous, piecewise linear function (on each interval Ij = [xj,Xj+i)) 
such that P\u h (xj) — Uj, and the piecewise constant interpolator is defined as 
Pou h (x) — uj for x E Ij. In view of the estimate ||D+u"||2 < a(t) and the fact that 

\\D + u% = \\P D + u% = \\d x P 1 u h \\ 2 , 

we see that there exists a function u <G Lj^ c ([0, oo); H 1 (R)) such that Piu h — > u in 
L\ oc {[Q,oo) x R) and d x P x u h A d x u in L°°([0, oo); L 2 (E)). Note that this proves 
the last assertion of the Lemma, since D + u h — d x P\u h . With this in mind, we 
wish to prove that u h = Pqu 11 — > u in L\ oc {R. x [0, oo)); For this, let fi be a bounded 
set of R x [0, oo). We have 

H^ou" - «IUi(n) < \\P^ h - Piu h \\ LHn) + \\P%u h - u\\ L i {n) . 

We have just seen that the second term above vanishes as h — > 0. For the first term, 
note that since O is compact, there are some M, Tel independent of h such that 

\\P u h - P lU h \\ LHn) < [ V f \P u h (x,t)-P lU h (x,t)\dxdt. 

JO i -i ,,,,, Jij 



|j|<M/h- 



Since 



y - P u h \dx = -\u j+1 - u^, 



we obtain from the estimate H-D+w' 1 !^ < a(t) in Proposition 3.4 and Young's in- 
equality that 

{■T 



\P Q u h - Piu h \\ L i{n) < f Yl 
|i|<M//i 



\j\<M/h 

dt 

|j|<Af/ft 362 

T 

< C / h(M + a(t)) dt 
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which tends to zero with h. This proves that Pqu — > u in L\ oc (K x [0,oo)) and 
completes the proof of Lemma |4.2| □ 



Proof of Theorem 2.2 It onl y remains to check that the (strong) limit (u,v) of 
(u h , v ) given by Lemma 4.2 is the unique entropy solution to the Cauchy problem 

( |l~2l^ -( |Oc| ). 

First, consider the Schrodinger equation (1.2a I. Taking the discrete equation 



(2. 7a I, multiplying by a test function 9 as in Definition 2.1 and integrating gives 
-iu h d t 6 + BD\u h - (\u h \ 2 u h d + g(v h )u h )6 dxalt = 0. 



// 

J JR 



x[0,oo) 



From the convergence properties in Lemma 4.2 we see that all terms except the one 
with D 2 u h converge to the corresponding ones in the equation (1.2a). Regarding 
that term, we easily find from the definition of the discrete derivatives 



// 

J jr 



x[0,oo) 



6D 2 h u h dxdt : 



D + u h 6' h dxdt, 



x[0,oo) 



where 6' h (x,t) :— [j-j. 1 0(x,t) dx — fj 9(x,t) dx)/h. Since 8' h converges strongly to 
d x 8, Lemma |4.2| ensures that this term also converges to the corresponding one in 
flTlal ). 

Next, we wish to prove that v is a solution to the conservation law (1.2b), in 
the sense of Definition 2.1 To this end, we must prove that for every non- negative 
smooth function 4> with compact support in K x [0,T), the inequality 



(4.17) 



V(v)d f ^ + ( gi (v) - q 2 (v)\u\ 2 )d x 

x[0,oo) 

(r 1 '(v)g'{v)-q 2 (v))d x (\u\ 2 )cj)dxdt+ [ V {v {x))<!>(x, 0) dx > 



holds. What we have is (cf. (4.3)) 
(4.18) 

r){v h )d t <t> + ( qi (v h ) - q 2 (v h )\u h \ 2 )d x 0dxdt 

x [0,oo ) 



and so 



(4.19) 



poo pOO 

- h^fyri'^dtVjdt- ^2<j) j+1 (q 1 {v j+1 ) - qi{v 3 ))dt 

I ^Z ( t ) 3+i{ c l2{v 3 +i)\u J+1 \ 2 - q 2 {v J )\u j \ 2 )dt- I 7](vo(x))(j>(x,0)dx, 
Jo -^n, Jr 



r)(v h )d t <t> + ( qi (v h ) - q 2 {v h )\u h \ 2 )d x (t>dxdt 

Rx[0,oo) 

A + B + C + D - I r](v$(x))(t>(x,0)dx, 



where A, B, C and D are defined before (4.4). 

From the strong convergence of {u h ,v n ) to (u,v) given in Lemma 
that the left-hand side of (4.19) converges to the first line of (4.17). Since (by 



4.2 



we see 
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assumption) Vq — > vo strongly, it it only necessary to check that 



(4.20) hm (A + B + C + D) > 



Rx[0,oo) 



4>{x,t){r 1 '{v)g'{v)-q 2 {v))d x \u\ 2 dxdt 



for (4.171 to hold. 



Consider the terms A + B. From the decomposition (4.5 1, the positivity property 



(4.6 ), the estimate (4.11 ) and the discussion after (4.11 ), we see that lim^o A+B > 
0. 



Let us now analyze the terms contained in C + D, see (4.12). First, we claim 
that lim/j^o Si > 0. Recalling (4.13), the property (4.14) shows that Su has a sign. 
So we must prove that 

22^+1 / rf'^duj+i^G^+ivj^j+i) - Gj,+) dvdt -+ 



Sl2 '- 

as h — > 0. We have 



Si2<c|M|oo I* E l^+i-^llK+il 2 -Kfl di 

J ° \i\<J 

< c ( f E %+i-^| 2 *) 1/2 ( f E rlK+il 2 



,1/2 



lil<J 



lil<J 

and, from ([3~9| , ( [37l3| ) , 

E ^|K + il 2 - M 2 | 2 < WAL E - ^i 2 ^ C ( T )- 



/l 1 

|3|<J 

Therefore, S3.2 — ^ if 



or, in particular, if 
(4.21) 



5>i 

lil<J 



Vj+l 



dt^O 



,h\ 2 



dxdt = \\v h -v h 1,2 



L 2 (0) 



where we have set v + (x) = v(x + h) and £1 = (0,T) xl. A standard result in 
compactness [131 P- 4] implies that if ||?;' 1 ||l 2 (£"2) = ll w +llL 2 (n)i an d if v h — > w in 
i 2 (57), then i/£ — ^ v actually gives — > v, and so (4.21) will hold. Therefore, we 
only have to prove the weak convergence vh. — * u in L 2 (f2). But if y € L 2 (S1) then 

ip(t,x)(v h (t,x + h) ~ v(t, x)) dxdt 

<p(t, x)(v h (t, x + h)- v h (t, x)) + tp(t, x)(v h [t, x) — v(t, x)) dx dt 
{ip(t, x — h) — ip(x))v h (t, x) + (p(t, x)(v h (t, x) — v(t, x)) dx dt 

<\\<p+- Mvih + yh\\v h - v\\ 2 -> 0. 



This proves (4.21 ) and so the term 5 12 goes to zero with h. 
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Next, we estimate the term S 2 given in (4.15|. We find 
S 2 = 



Oi) K, 2 >(,(</, • G j>+ ) 



+ {<t>j+i - <f>j)"nj(\uj\ 2 Gj,+ - Gj,+) - <f>j+Wj9j{\uj+i\ 2 - W 3 \ 2 ) dt 



S- 



21 T O22 T" 023- 



Thus 



S21 < /ill^'IU / Y] W 3 \ 2 Vj{ - 9j - Gj.r 



\j\<J 



and so the same computation as the one after (4.15) applies, giving 

S21 -> 0. 



Next, using the estimates (3.7) and (3.9), 



S 2 2<Ch\\cf> , \\ OD f Y,W u j+i\ 2 ~ Kf|* 
Jo 



< c(T)h\\4>' w^^o 



as h — >■ 0. Finally, from Lemma 4.2 and from D + \u \ — v c^up (which is a 



consequence of the same lemma), we see that 



'23 



<f>(x,t)T]' \v)g' \v)d x \u\ 2 dx dt. 



It only remains to analyze the terms D 3 and C 2 in (4.12). We have, again from 
Lemma 14.21 



/>oo 

D 3 = ^4> ]+1 q2{v 3 )(\u J + 1 \ 2 - 

JO Sr .~ 



uA 2 ) dt 



4>(x, t)q2(v)d x \u\ dx dt 



as h — » 0. For the term C2, we easily reduce it to the cases already treated by 



proceeding as in the comments after (4.16), giving C2 — > 0. This completes the 
proof of Theorem 2.2 □ 



5. Numerical experiments 
In this section, we present some numerical computations illustrating our results. 



5.1. A fully discrete scheme. To implement numerically the scheme (2.7), we 



choose a semi- implicit Crank-Nicholson scheme for the Schrodinger equation (2.7a), 



with a Newton iteration for the nonlinear term |u| 2 w, coupled to a semi-implicit 
Lax-Friedrichs scheme (see (|2.8|)) for the conservation law (2.7b). More precisely, 



24 



P. AMORIM AND M. FIGUEIRA 



given a spatial mesh size h and a time step r, we consider the following algorithm: 
(5.1) 

- + «?)| 9 i(«y +1 +«?) + «?), 

+ ^'K + i)K + il 2 - ^(«7_i)l«?-il a ) + ^Ki 1 - K +1 + "ft 1 ) 

+ 2^[ ( ^ - W + |U ^ 1|2) - - ^ + i)^(l^| 2 + l«?-il a )]- 

Of course this scheme, being at most first order, serves only to illustrate our results. 
Indeed, it is well known that the Lax-Friedrichs scheme suffers from a degradation 
of solutions due to numerical viscosity. More accurate discretizations would be 
interesting to implement, for instance, by considering more powerful time-stepping 
methods, or by implicitly coupling the two equations. 

As is usual in simulations involving the Schrodinger equation on the whole line, it 
is reasonable to consider the problem in a sufficiently large domain, and taking care 
that the value of the computed solution remains very small (in our case, this value 
is usually of the order 10 -7 after a few thousand iterations) near the boundary. 

In all simulation presented below, the computed L 2 norm of the solution u h 
is exactly conserved (up to a very small error committed in its computation), in 



accordance with the L norm conservation property (3.2) 



5.2. A test case with linear flux. As a first example, and to test the accuracy 



of the numerical scheme (5.1), we consider the following problem in which the 
conse 

(5.2) 



conservation law has a linear flux, 

{idtu + d xx u = k(vu + g|w| 2 u) 
d t v + 1 d x v = sd x (\u\ 2 ) 

whose exact traveling wave solution can be found (see, for instance, 




u(x, t) = e tM e lc{x ~ ct)/2 J ~ sech (Ve(x - ct)) , 



, ^ 2E , 2 

v{x, t) = a— —sech 



(5.3) 

d(t. f\ = n, 

l/3| 

where c, A > are to be chosen, E := A — c 2 /4, a := 8/ (7 — c) and (3 = k(a + q). We 
must have E > 0, j3 < to ensure existence of a solution (see [16] for the details). 
Here, we choose q = k = X = j = l, a = — 2 and c = 3/2. 

Note that here and in the following examples, we have omitted the function g(v) 
from the formulation of the problem. This is because we are taking g'(v) to be 
the characteristic function of some large interval [—M,M]. In all our numerical 
experiments, the value of \v\ never gets close to M, and thus among the tested 
initial data, the original problem can be stated without reference to g. 

In Figure [l] we present the initial data \uq\ and Vq corresponding to (5.3) with 



t = 0, together with the computed solution \u h \, v h , which can be observed to be a 
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Linear flux case, 6000 pts, dx= 0.01 , T = 5 
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Schrodinger Initial data 
Conservation law Initial data 
Schrodinger solution 
Conservation law solution 



Figure 1. Initial data and computed solutions of system (5.2 1 



Square of L2 error, linear flux, dt= 0.01 , T = 3 
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Figure 2. L 2 error between the exact solution (5.3 1 and the com- 
puted solution, linear flux. 



traveling wave. The error in the L? norm between the exact solution ( |5.3[ ) and the 
computed solution for various spatial discretizations is presented in Figure [2] 



5.3. A test case with nonlinear flux. We now test the scheme ( |5.I[ ) on an 

explicit solution to the problem with nonlinear flux f(v) = v 2 but in which the 
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Square of L2 error, Nonlinear flux, dt= 0.01 , T = 3 
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Figure 3. L 2 error between the exact solution and the computed 
solution, nonlinear flux. 



first equation is linear in u (we thank J. P. Dias for bringing this solution to our 
attention): 

{idt.u + d xx u = vu 
d t v + d x v 2 = d x (\u\ 2 ). 

The explicit solutions can be found by the following Ansatz, 

u(x,t) = e lbt r(x), 
v(x t t)=p{x). 

From the equations we get 

/3 = —\r\ and - br + r" = -r 2 , 

of which a solution is known to be r(x) = 6(3/2) sech 2 (v / &x/2). Thus 

(u,v) = (e M r(x),-r(x)). 

In Figure [3] we present the error in the L 2 norm between this exact solution and 
the computed solution for various spatial discretizations. 

5.4. The general case. Here we present some numerical solutions to the full sys- 
tem ( 1.2a I ,( 1.2b I , with f(v) — 3v 2 , showing (the modulus of) a solution of the 
Schrodinger equation interacting with a solution of the conservation law. 
We take as initial data the functions 

Uo(x) = e 5lx/2 V6sech{V3x), 

v (x) = X[-10,10] 

on the spatial domain [—50, 50]. In Figures [3jj6] we can observe the usual behavior 
of the solutions to a nonlinear conservation law, such as the propagation of the 
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initial shock discontinuity and the formation of rarefaction waves. In addition, the 
interaction between the two equations induces the formation of new waves, as can 
be seen clearly in Figure [5] 

For completeness, we show in Figure [7] the real and imaginary parts of the 
solution u of the Schrodinger equation, at time T = 2.5. 

Again, in these simulations, there was no need to consider the coupling function 
g{v), which is set to a cutoff function for some large cutoff parameter, since the 
value of the solution v never approaches this cutoff value. 

Recall that one of the purposes of the coupling function g was to ensure an L°° 
bound on the solution v, which is in turn essential to our convergence (and well- 
posedness) proof. Such an L°° bound is an open question when g(v) = v [21 [9]. 
Now, in our experiments, we observed no blow-up phenomena in v. Since in our 
simulations the solution u was always very regular, this suggests that, when the 
function g is the identity, any eventual blow-up of v should be associated with 
initial data for u (and thus source terms for the conservation law (1.2b)) having 
the minimum allowed regularity, namely u G H 1 (M.) but u x L°°(M). A thorough 
study of this case, however, is out of reach of our simple code, and would need a 
more sophisticated numerical approach. 
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Figure 4. Initial data and computed solutions of system ( 1.2a), ( 1.2b ) 
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Full system, 5000 pts, dx= 0.001 , T = 1 
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Figure 5. Computed solutions of system ( 1.2a), ( 1.2b I, T = 1, 
T = 1.5. 



[17] F. Murat, Compacite par compensation, Ann. Scuola Norm. Sup. Pisa CI. Sci. (4) 5 (1978), 
no. 3, 489-507 

[18] L. Tartar, Compensated compactness and applications to partial differential equations, 
inNonlinear Analysis and Mechanics: Heriot-Watt Symposium, vol. IV, pp. 136-212, Pit- 
man, Boston (1979) 

[19] E. Tadmor, Numerical viscosity and the entropy condition for conservative difference schemes, 

Math. Comp. 43 (1984), no. 168, 369-381 
[20] M. Tsutsumi and S. Hatano, Well-posedness of the Cauchy problem for the long wave-short 

wave resonance equations, Nonlinear Anal. 22 (1994), no. 2, 155—171 



P. AMORIM AND M. FIGUEIRA 



Full system, 5000 pts, dx= 0.001 , T = 2 
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Figure 6. Computed solutions of system ( 1.2a), ( 1.2b I, T = 2, 
T = 2.5. 
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Full system, 5000 pts, dx= 0.001 , T = 2.5 
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Figure 7. Computed solution of system ( 1.2a), ( 1.2b ). 



